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The correlated Kondo-lattice model is used to describe the interaction of electrons in a single 
conduction band with localized magnetic moments as well as their mutual repulsion. It is our 
intension to provide an analytical exact result for this model by considering a two-site cluster. 
An equation-of-motion approach as well as the spectral theorem is chosen to obtain the complete 
expression for the one-particle Green's function. In a forgoing article the problem has been solved 
for an insulator, the present paper is devoted to the technical more demanding task of a metal. 



I. INTRODUCTION 

The magnetism of a large variety of so-called local mo- 
ment systems, like manganites, gadolinium or europium- 
chalcogenides, is governed by an interplay of the proper- 
ties of magnetic moments localized at certain atoms and 
itinerant conduction electrons propagating through the 
lattice. If the system of localized spins Si were treated in- 
dependently, a certain magnetic order could be explained 
with the help of a Heisenberg-like exchange coupling 
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In contrast to that the Kondo-lattice model tries to ex- 
plain magnetic order via an indirect coupling of Si and 
Sj mediated by the spin cr of the itinerant conduction 
electrons. The importance of the on-site interaction 
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was first emphasized by Zener0. 

As far as the conduction electrons are concerned they 
might be described by their kinetic energy 
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only. Here, c, ff is a fermionic annihilation operator in 
second quantization. However, several authorsa pointed 
out that at least for manganites with its relatively narrow 
e s -bands an ansatz that neglects any electron-electron 
interaction is an insufficient treatment. This is because 
U is probably the largest energy scale in the problem. 

Another indication of this insufficiency is the exact 
solution of the zero-bandwidth limit (Ty = T 5ij, also 
called atomic limit) of the Kondo-lattice model H s +H s f- 
In this limit the one-particle Green's function has an ex- 
citation spectrum of four energy poles 
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where S is denned by S? = h 2 S(S + 1). If the "test" 
electron enters an empty lattice site the energies e\ or 
£2 are necessary. The other two energies occur if the 
entered lattice site is already occupied. Apparently £3 
has the lowest energy for J > 0. Because of the involved 
Coulomb repulsion this is not understandable. 

The simplest way to include electron-electron interac- 
tion in the Hamiltonian is an on-site Hubbard interaction 
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It has been showni that this has an significant effect on 
the spectral function of the Kondo-lattice model. For 
the atomic limit calculation it implies that the energies 
£3 and £4 are shifted by U and £1 becomes lowest. 

The combination of all these possible effects in one 
Hamiltonian 



H = H s + H s 
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is in this paper called correlated Kondo lattice model 

(CKLM). Due to the fact that there is no analytical solu- 
tion of the CKLM available, it is essential to obtain exact 
statements from liniiting,cases. 

In a foregoing articloa we explained, why especially 
the generalization of the atomic limit result to a two-site 
cluster is important for a better understanding of the 
CKLM. There, we investigated the Hamiltonian H s +H s f 
with an assumed hopping integral of the form 

!T i j and i, j in same cluster, 
T i = j, 
i,j in different clusters. 

The part Hff was omitted to gain the chance for a con- 
nection of the cluster to the lattice. This is possible if 
expectation values of spin operators, like \Sf) are not 
uniquely determined by the properties of the cluster. 
Furthermore, a Hamiltonian without Hff is much eas- 
ier to handle. The part H ss , on the other hand, was 
redundant since the limit of an empty conduction band 
was considered. For this situation, which is related to 
an insulator, we were able to give the complete analyt- 
ical expression for the one-particle Green's function as 
obtained by a solution of the equations of motion. 
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With the present paper we try to obtain similar results 
for metals. Since the occupation of the conduction band 
is now arbitrary, the mathematical systems that need to 
be solved are about one order of magnitude larger than 
in the case of an insulator. Additionally, it is necessary 
to take the Hubbard interaction into account. Therefore, 
we consider the Hamiltonian 
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where 5 T = S+,S i = S~ and z T = +1. zj_ = -1. For 
simplicity we restrict ourself to 5 = 1/2. 

Despite of the larger effort, we use the same techniques 
already proposed for the insulator. What has been ex- 
plained on the more concise problem there, can now be 
applied directly. It will not be possible to provide ev- 
ery detail of the calculations here. A thorough docu- 
mentation, of all necessary transformations can be found 
elsewhereB. 



II. SOLUTION FOR THE EQUATIONS OF 
MOTION 



To describe the properties of the electronic partial sys- 
tem of the cluster we are especially interested in the one- 
electron retarded Green's function ((c i(r ; c]- CT )) B . This 
function is closely related to the spectral density and 
contains all information about possible single-particle ex- 
citations. 

The method of choice is an equation-of-motion tech- 
nique. Compared to a direct determination of the 
Lehmann representation it has two advantages. First, 
an equation of motion 
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contains as an inhomogeneity a correlation function 
([A,B} + ). This will directly appear in the numerator 
of the Green's function, without being forced to actu- 
ally calculate the expectation value. Therefore, we are 
able to discuss the dependence of the spectral weights 
on parameters as the averaged band occupation or the 
averaged z-component of the localized spins. Secondly, 
each equation of motion contains higher Green's func- 
tions (([A, H]-; Bfj E . Hence, by using intelligent tech- 
niques we do not only get ((A; B^j E but a complete set of 
Green's function with much more information. We will 
make use of this fact in the next section. 

Since we are dealing with a finite size system, the hier- 
archy of equations of motion necessarily has to terminate. 
The complete set of 1040 equations can be formulated as 
a single matrix equation. By unitary transformations 



this matrix needs to be reduced to blocks of a size which 
can be solved directly. A first step in this direction is 
the implementation of the site symmetry. Therefore, we 
define combined Green's functions 
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where the site index s represents the opposite of site s 
(s=l=>s=2;s = 2=>s = l), C s and D g are arbitrary 
products of spin- and Fermi-operators. The two partial 
systems for e = +1 and e = — 1 can be solved separately 
and simultaneously (they only differ in the sign of T) . 

Even if additionally the mirror symmetry with respect 
to the spin quantization axis is used, one is still left with 
the large number of Green's function listed below: 
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Here, 5 stands for all possible products of spin operators. 
If, as assumed, 5=1/2 there are 10 such possibilities. 
The right column gives for each Green's function the cor- 
responding correlation function emerging in its equation 
of motion. It is useful to group these functions into four 
density classes, determined by the number of fermionic 
creation operators in its active part: 
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The Green's functions of density class IV vanish whenever 
the number of participating electrons is smaller than 3. 
The ones in density class 11 vanish if we have less than 
2 electrons etc.. Using this notation the matrix equation 
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of all Green's functions has the block structure 
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where the right hand side represents the column of all 
inhomogeneities appearing in every equation of motion. 

Since the matrix has a triangular shape, the remaining 
task is to find the inverse of the blocks Mi, Mj, Mm and 
Mw- K is therefore possible to stay in certain density 
classes, which justifies the grouping. Additionally, we 
can benefit from particle-hole symmetry. It implies that 
the density classes I and IV as well as the classes I and 
H can be treated in a similar way, respectively. 

The density class I consists of the 10 Green's functions 
with only one fermionic operator but with several spin 
operators. If one could assume that all Green's functions 
of the other density classes vanish, one just had to in- 
vert the matrix Mi to get a solution. The assumption 
is fulfilled if and only if one describes an insulator with 
a single excess electron in a-ptherwise empty conduction 
band. The foregoing articlea was devoted to these mate- 
rials. In the case of arbitrary band occupations we have 
to solve the equation Mi ■ G\ = K\ — M\$Gy.. Since 
only the right hand side is modified we can use the same 
rules for the combination of spin operators in the Green's 
functions as in the case of an insulator to simplify Mi. 

Matrix Afn is also simplified by an intelligent combina- 
tion of Green's functions. As far as the spin operators are 
concerned we use again the rules proposed in the treat- 
ment of an insulator. To get an idea about the behavior 
of the fermionic operators we investigated another lim- 
iting case: If in H the parameter J is set to we have 
a system without any spin operator, namely a Hubbard 
cluster. Its solution is sketched in appendix [A| The rules 
for a formation of combined Green's functions used there 
can then be applied to equation ( ^f]| ) . 

This technique to combine Green's functions by a suc- 
cessive treatment of the spin operators and the fermionic 
operators respectively works, but is rather lengthy and 
fault-prone. When doing this job we learned that 
a simple text editor can be awesomely helpful. Al- 
ready when generating the 1040 equations of motion a 
cut-and-paste mechanism can be used effectively, since 
the same commutators [c sa ,H]-, [5f,W]_ and [Sg,H]- 
have to be implemented several times. For combinations 
the search-and-replace function is quite helpful. 

Appendix [B| gives detailed instructions which Green's 
functions need to be summed and subtracted to achieve 
a simpler structure of Mn. If implemented one obtains 
matrix blocks of size 3x3. There are four types of these 
matrices: 
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Here, we used the abbreviation a = | j and [i\ 
are sign parameters. Mo, which already played a role in 
the insulator problcmQ, has the eigenvalues 
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The other matrices are connected to non-vanishing 
electron-densities as can be seen at the appearance of 
U. The structure of Mi is simple, and its eigenvalues are 

E[^=^U, E 12 = 2a, E 13 = -2a. (15) 
In contrast to that the determinantal polynomial of M 2 

= E 3 + fi 5 E 2 U - AE (a 2 + T 2 ) - 4a 2 /i 5 [7 (16) 

can only be factorized with Cardan's formulas. The three 
eigenvalues of this cubic equation are 
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Similarly, we have for M3 the cubic equation 
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where C% = 
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If the eigenvalue problem of the cluster is studied and 
especially the Hilbert subspace of two electrons in the 
conduction band is considered, one is confronted with ex- 
actly the same two cubic equations ( p"6] ) and (j2C|). This 
underlines the close relationship of the eigenvalue prob- 
lem and the equation-of-motion technique. 

To determine the inverse of M x , x = 0, . . . , 3 it is suf- 
ficient to know its eigenvectors. As explained more de- 
tailed in the treatment of the insulatorta, basic algebra 
then immediately allows to give an expression for M ~ : 
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with well defined coefficients m xk (i, j). This structure is 
particularly suitable because we intend to get the ana- 
lytic expression for the one-particle Green's function as 
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with p linear energy poles. Since we succeeded to find 
combined Green's functions which are connected in 3 x 3 
systems, an inversion of the matrices like in (|^) will en- 
sure that every combined Green's function has the form 
(p5|). To return to the original Green's functions one just 
has to undo the combinations. That happens by summa- 
tions and subtractions of the combined expressions and 
will therefore maintain the form (p5|). 

We mentioned above that the calculations can be done 
within certain density classes. To be more precise, to 
solve equation (10) the following four matrix equations 



have to be considered: 
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One has to obey this order, because apart from system 
IV we always have an input of a class with higher electron 
densities. On the one hand, the matrix multiplication on 
the right side implies an additional summation of already 
calculated Green's functions. On the other hand, a par- 
tial fraction expansion is necessary to return to the form 
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Occasionally, it happens that Em = Em, leading to 
quadratic energy poles. However, after performing all 
necessary summations the spectral weight of these poles 
will always vanish. 

Taking all this into account, one obtains after a 
straightforward calculation expressions for the Green's 
functions of density class I. If furthermore all combi- 
nations are undone, the desired analytical exact result 
for the one-electron Green's function Uc s • f> 
given. It consists of 102 energy poles. 

Six of them, namely the energies 
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are a result of the inversion of Mi. They are just a subset 
of the energies 
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which occur, if Mn is inverted. The structure of these en- 
ergy poles is quite clear. The Green's functions of density 
class I include excitation processes from a one-electron 
state to a two-electron state. The energies of the first are 
mediated by the eigenvalues of Mo, the energies of the 
latter are connected to the matrices Mi , M^ and A/3 . To 
is just a trivial energy shift, —a actually belongs into the 
definition of Mo. The super- index "±" in Eo k illustrates 
that always both possible signs of T need to be consid- 
ered, as mentioned subsequently to (0). Hence, Mj leads 
to 3x3x3x2 = 54 energy poles. 

It has already been pointed out, that Mm can be 
treated similarly to Mj according to the particle-hole 
symmetry of the system. This symmetry can be inves- 
tigated more thoroughly when the Hamiltonian is refor- 
mulated in terms of hole-creation and -annihilation op- 
J _ = z„C;„. It leads to the result that with 
f(T, J, U) of the one-particle 
To - f{-f, J, U) is an energy 
2|) one gets all the energies 
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we have six degenerate energies, leading to p = 2 x 54 — 
6 = 102 energy poles of the one-particle Green's function 
((c i(T ; c| CT ))„ and all higher Green's functions. 

All these excitation energies are functions of the model 
parameter T, U and J. Since the cluster is a generaliza- 
tion of the atomic limit, the T-dependence is of special 
interest (see figure [l]). Because of the large number of 
energy lines, the formation of two Hubbard bands can 
clearly be seen. Although they are well separated at 
small T, there seems to be an overlap of the bands above 
0.5 eV. When going to higher values of T a separation of 
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Hopping-Integral T [eV] sf-coupling constant J [eV] 



FIG. 1: Dependence of every energy poles of the one-particle 
Green's function on the hopping integral T. The fixed model 
parameters are: To = OeV, U = 2eV and J = 0.2 eV. 




Hopping-Integral T [eVJ 



Hopping-Integral T [eVJ 



FIG. 2: Dependence of the energy poles on the hopping in- 
tegral T. At the left the behaviour of the CKLM is shown, 
for comparison the Hubbard cluster (J = 0) is investigated at 
the right for the same parameter range. 



the bands and a formation of bundles occurs. As can be 
seen when the two parts of figure |^ are compared each 
of these bundles results from the J-splitting of an energy 
pole obtained for the two-site Hubbard cluster (appendix 
[A|) . This splitting can best be shown, if, as in figure [|, 
the J-dependence of the energy poles is studied. The 
picture suggests a symmetry with respect to the sign of 
J, but this is not completely given. 

Coming back to figure ^ another striking fact needs 
to be mentioned. It concerns the limit T — > 0, which is 
identically to the situation of the atomic limit. As noted 
in (|J) we expect in this limit a reduction to four energy 
poles. In contrast to that we obtain twice as many energy 
lines for T = in figure [l]. Obviously, the energies alone 
are insufficient to describe excitation processes. Remarks 
on the corresponding spectral weights are necessary. 

The routine we suggested for the solution of the equa- 
tions of motion inevitably gives for each Green's function 
these spectral weights, too. The fact that we don't give 
the complete analytical solution here, has two main rea- 
sons. First of all, the routine consists of quite a lot of 



FIG. 3: Dependence of the energy poles on the sf-coupling 
constant J. The fixed model parameters are: To = OeV, U = 
2eV and T = 0.4 eV. The vertical line marks the J value of 
figure [[]. 



complicated summations: The combined Green's func- 
tion for instance are just sums. At least one summation 
is necessary to get the inverse of a single matrix as given 
by equation (24). Their multiplication with a column of 
correlation functions as well as the matrix multiplications 
in each of the lines of (|2^) imply further sums. At the end 
partial fraction expansions have to be performed and the 
combinations of Green's functions needs to be undone. 
Even so, all this is simple arithmetics it leads to really 
lengthy expressions for the spectral weights. Secondly, 
the spectral weights of each Green's function depend in 
principle on all correlation functions that appear in in- 
homogeneities of the same or higher density classes. For 
((c^; c \i)) E this means a dependence on 1040 correlation 
functions (1/3 is zero by definition). All of them need to 
be specified. Only afterwards the result for the spectral 
weights is really meaningful. 



III. TREATMENT OF THE CORRELATION 
FUNCTIONS 

Already when dealing with the limiting case of an 
empty conduction banco one has to develop concepts for 
the determination of the emerging spin correlation func- 
tions. On the one hand, a simpler Hamiltonian, that 
describes the spin system and includes lattice effects like 
anisotropy or molecular fields, can be investigated. Al- 
ternatively, the spectral theorem provides information on 
these correlation functions. The latter concept implies 
that all calculations are based on the cluster Hamilto- 
nian only and no connection to the lattice is made. In 
this paper we are going to use a combination of these two 
concepts. 

As mentioned before, it is a crucial advantage of our 
methods that besides ((c ia ; c \ a )) a complete set of all 
Green's functions participating in the equations of mo- 
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tion is obtained. For an arbitrary Green's function, by 
construction given in the form fl25|), the spectral theo- 
rem allows the calculation of a corresponding correlation 
function 
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where (3 = (fcsT 1 ) - . Therefore, we obtain expressions 
for 1040 correlation functions (again 1/3 is zero by def- 
inition). Their principle structure is equal for Green's 
functions of the same density class. Table | allows a 
comparison with the structure of the corresponding in- 
homogeneities. 
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TABLE I: For Green's functions of every density class the 
principle structure of the correlation functions emerging in the 
equation of motion (B) and obtained by the spectral theorem 
(K4j) is given. 

The application of the spectral theorem to Green's 
functions of density class IV results in correlation func- 
tions that do not appear in the spectral weights we are 
looking for. Only the application to the 960 Green's func- 
tions of the remaining density classes is helpful. For each 
of them the spectral theorem provides an equation that is 
linear in the correlation functions, because this linearity 
is given for the spectral weights a k in (|34|). In principle 
we can therefore write a homogeneous matrix equation 
for the column of correlation functions Ky. 
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Written in this form, equation ( |35| ) is underdeter- 
mined. This seems to offer the possibility to use addi- 
tionally the other concept for the correlation functions, 
that is to make assumptions for certain expectation val- 
ues or to specify them on the basis of a Hamiltonian that 
includes lattice effects. To achieve a non- vanishing mag- 
netic order of the system of localized spins it is necessary 
to specify the spin correlation functions, {S). Lattice 
effects can be included if these expectation values are 
based on the Hamiltonian Tiff +Hf with 



H, = -b(Si + r,SS), 
as discussed in the forgoing article 
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Whenever a certain correlation function is specified, 
the matrix equation ( |35| ) becomes inhomogeneous: 
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With the inhomogeneity in the equation of motion of 
(( c io-! c L))b' which is 1 by definition, at least one speci- 
fication is inevitable, enforcing a nontrivial result for the 
set {Ki},. 

However, the freedom for the specification of certain 
correlation functions is drastically reduced because of 
additional information available. We have already men- 
tioned that 1/3 of the K^s is zero by definition, which 
reduces the number of unknown variables in ([55]). Even 
so, in 1/3 of the cases the application of the spectral 
theorem (|4]) does also give zero at the left hand side, 
this does not lead to a reduction of the number of equa- 
tions since the right hand side is still a linear relation of 
the variables. Furthermore, we know that certain pairs 
of correlation functions have to be identical due to sym- 
metrical reasons. 

Hence, we actually do not have an under-determined 
problem, but rather more equations than variables. The 
indicator for a proper choice of equations is a nonzero 
result for the determinant of the matrix in (|3?j). Rules 
can again best been studied, if the limiting case J = 0, 
the Hubbard cluster, is considered. 

For an overdetermined problem the introduction of 
additional specifications for correlation functions is of 
course problematic. They are likely to contradict some 
other information we have about the system. For the sake 
of a connection between cluster and lattice we want to ac- 
cept that. The consequences shall be demonstrated with 
an example: If the atomic limit (T = 0) is considered, 
7 nontrivial correlation functions appear in the spectral 
weights of the one-particle Green's function. Since the 
spectral theorem provides only 6 equations, (Sf ) can be 
specified on the basis of a Brillouin function for the lat- 
tice. However, if a non- vanishing value for (Sf) is cho- 
sen one has to be aware of the fact that the equality 
(S'J t "cJ|C 1 |) = (5 , ]/ c iT c ii) w * u be broken. The equality 
itself follows from the assumption that expectation val- 
ues are real, but can also be derived from particle-hole 
symmetry. 



If an equation like (37) needs to be solved exactly for 
the correlation functions {Ki} { standard algebraic tech- 
niques like Gauss' elimination method can be used. In 
our work we implemented all these tedious manipulations 
in a computer program. Its input values are the model 
parameters and assumed specifications for the correla- 
tion functions. By setting successively one and only one 
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correlation function to 1.0, calculating the set of Green's 
functions (here the analytical exact expressions are used) 
and applying the spectral theorem, the matrix ([35]) is de- 
termined. Afterwards the matrix equation (|37| ) is solved 
numerically exact. 

The last step contains a specific numerical problem. 
Even if a proper set of equations is selected, for large 
ranges for the chemical potential /i the determinant of the 
matrix in ( |37| ) might be rather small. This is connected 
to the fact that the Fermi distribution function in (|4|) is 
almost a step function for low temperatures. The effect 
can already be studied when looking at the atomic limit 
for the Hubbard model 



ft = T h c + 7T K^-a- 



(38) 




c.p. [eV] 



E [eV] 



After the application of (|34| ) the two equations (a =f, |) 

(h a ) = (1 - (n_ a )) /_(T ) + (n_ a )f-(T + U) (39) 

define the correlation functions. Whenever n is located 
well between To and Tq + U and additionally k-sT <C U, 
it is /_(T ) - /_(T + U) « 1. Then (p|) becomes 




f-(T 




(40) 



with a vanishing determinant of the matrix. Ana- 
lytically, one can strictly show that (n*j = = 

/_(To)[l + /-(To)-/_(T + C/)]~ 1 . The numerical 
problem can be avoided if the chemical potential is close 
to an energy pole or if the temperature is high enough. 

For the two-site Hubbard cluster Ti s + H ss the 
mathematical systems are much smaller. Hence, al- 
ready temperatures around 100 K are sufficient to 
get meaningful results (see figure [I| ) over the whole 
range of the chemical potential (c.p.). In regions for 
/j, which are not close to the energy poles all spec- 
tral weights remain constant. Figure ^ also nicely 
illustrates a further consequence of the particle-hole 
symmetry. As already mentioned above, the energy 
poles of the one-particle Green's function occur in pairs 

(E k = T + f k (f, J, U), E k = To + U - f k (-f, J, 17)) . 

We can now see that for the corresponding spectral 
weights the implication 

& k = ot k {^,a,f) =>• a k = a k (U - /j,, -a,-f), (41) 

that can be proven analytically, is fulfilled. The signs 
of T and a are not essential because of site and mirror 
symmetry, respectively. Instead, the /z-dependence is the 
most important feature. 

Coming back to the complete CKLM we are again es- 
pecially interested in the dependence of the cluster prop- 
erties on T. For T_.= we are able to reproduce the 
atomic limit resultsu completely. This includes state- 
ments on the spectral weights and the fact that it van- 
ishes for the four additionally emerged energy poles in fig- 
ure [lj. The distribution of the spectral weight for higher 



FIG. 4: For the two-site Hubbard-cluster the dependence of 
the spectral weights of ((c iCT ; c| CT \^ on the chemical potential 
(c.p.) is given. The energy poles are situated at energies E 
given in appendix [X]. The model parameters are: To = J = 
0eV,f = 0.4 eV and U = 2eV. A temperature of T = 100 K 
is used. 




ll.'ppnie liueer.il T cV ] 



Energy poles E^ [eV] 



Energy poles 1^ [eV] 



Energy poles E^ [eV] 



FIG. 5: top left: as figure ^. others: Distribution of the spec- 
tral weight of (( c io-) C L)) B f° r different values of the hopping 
integral T (as given). The chemical potential (doted line) is 
always chosen such, that the determinant in ( p7| ) is not too 
small. All spin-correlation- functions are assumed to be 
zero. The other model parameters are: To = 0.0 eV, J = 
0.2 eV, U = 2.0 eV. Temperature: T = 500 K. 



values of T is shown in figure |j| Especially for T values 
higher than 0.5 eV the dominance of the energy poles of 
the Hubbard cluster can clearly be seen. 

The splitting of theses poles and a greater width of 
the bundles in figure ^ is connected to an increase of J. 
Nevertheless, a redistribution of spectral weight happens 
just gradually, as can be seen in figure [| It is remarkable 
that up to \ J\ = 0.25 eV the spectral weight is restricted 
to a few energy poles. One can also see that some of the 
energy poles get small negative spectral weights. This is 
caused by the two problems mentioned above. On the one 
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sf- 


oupling constant J [cV] 


J = +0.1 eV 









Energy poles E (c\ 



Energy poles E t [cV] 



Energy poles EJcV] 



J = -0.2 eV 

..li 1 


J„. . ,. .ii. 







Energy poles E [cV] 



Energy poles E (cV] 



FIG. 6: top left: as figure |^. others: Distribution of the 
spectral weight of ((c^; c ]f)) B f° r different values of the sf- 
coupling constant J (as given) with opposite signs. The chem- 
ical potential is kept fixed at /i = 0.0 eV. For the spin expec- 
tation values the specification (Sf\ = (Sf) = 0-1; (SIS^) = 
(S^S^) = (SiS?) = 0.0 has been assumed. The other 
model parameters are: T = 0.0eV,T = 0.25 eV, U = 2.0 eV. 
Temperature: T = 500 K. 



hand we have specified the values of the spin-correlation 
functions in figure ||. The choice does certainly not match 
with calculations based on the two-site Hamiltonian (|7). 
We have stressed above that this might lead to a break- 
down of certain symmetries of the cluster. The fact that 
every spectral weight must be positive, as can be shown 
with the help of Lehmann's representation, is a property 
that follows from such symmetries. On the other hand 
we also mentioned that an unfavourable choice for the 
chemical potential [i implies large numerical errors, due 
to very small determinants of matrices that need to be 
inverted. This can also cause negative values for spectral 
weights. 

Our program allows the presentation of the excitation 
spectrum for various situations. The main intension for 
the improvement of the atomic limit calculations was the 
ability to treat also an antiferromagnet order of the sys- 
tem of localized spins. An example for typical results is 
given in figure [?]. As in most of the other figures one 
can see, that even though the Green's function possesses 
more then 100 energy poles, only a few of them have a 
non- vanishing spectral weight. The expectation of quasi- 
continuous energy bands, provoked by pictures for the 
energy poles as figure [I], is not met by the result ob- 
tained from the calculation of the spectral weights. Like 
in the atomic limit we end up with a limited number of 



„ 0.2 - 

bis 

'53 

•3 O.lh 



-0.1 









' 1 

- 






J. .1 


- 

l v. J A 






1 T I 1 



Energy poles E [eV] 



FIG. 7: Distribution of the spectral weight over the en- 
ergy poles of (tyC^\c\-\)) E - By the specification of the spin- 
correlation functions (Sf) = 0.1, (Sf) = -0.1; (SfS%) = 
(SfS'^} = (S^S^) = 0.0 an antiferromagnetic order is 
assumed. The model parameters are: To = 0.0eV,T = 
1.0 eV, J = 0.2eV,£/ = 4.0 eV. Temperature: T = 500 K. 



delta peaks in the spectrum. 



IV. SUMMARY 

We considered the analytical solution for a two-site 
cluster described by the CKLM. In contrast to a foregoing 
publicationcl we did not restrict the occupation of the 
single conduction band. Moreover, a correlation effect in 
the form of a Hubbard interaction in this band is taken 
into account. The only remaining constraint is S = 1/2 
for the localized spins. 

The mathematical effort to achieve the analytical exact 
result for the Green's function ((c icr ; cj^))^ is immense. 
It was the goal of the paper to show that it is still man- 
ageable and to explain the necessary techniques. Major 
points were the separation of the 1040 Green's functions 
in four density classes, the use of intrinsical symmetries 
and the study of limiting cases as the empty conduction 
band and the Hubbard cluster. As a consequence we 
found the "correct" combinations of Green's functions 
(appendix |b|) to obtain system sizes of order 3x3. 

The suggested routine for an analytical calculation 
consists of two main parts. First of all the set of equa- 
tions of motion needs to be solved by: Combining their 
inhomogeneities correctly, including the results for higher 
density classes as in equation ( 26 ) , inverting the 3x3 ma- 



trices, performing partial fraction expansions and undo- 
ing the combinations of the Green's functions. Secondly, 
the correlation functions in the spectral weights need to 
be determined. This happens by: Specifying a desired 
subset of the correlation functions, applying the spectral 
theorem to every determined Green's function, selecting 
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a representative set of these linear equations, using an 
LU-decomposition to solve the linear system of equations 
for the unspecified correlation functions. 

The first part has been done by hand and the ana- 
lytical expressions for the obtained 102 energy poles as 
well as a graphical presentation of the excitation spec- 
trum has been given in this paper. The correlation effect, 
which leads to a formation of two Hubbard bands, and 
the characteristic dependence on the important model 
parameters can clearly be seen in figure || and |[ For 
the second part a computer program has been used, since 
there exist reliable routines for the numerically exact de- 
composition of large systems. The resulting vector for the 
correlation functions was used to discuss the distribution 
of the spectral weight among the excitation energies. A 
major conclusion from the provided figures is the concen- 
tration of the spectral on a few energy poles. The notion 
of almost "continuous" bands formed by 102 energies is 
not confirmed if the numerators of the Green's functions 
are considered. 

We propose to use a subset of the correlation functions 
in the numerators of the Green's functions to construct a 
connection of the cluster to the lattice. This freedom in 



the specification of expectation values strongly influences 
the results for the spectral weights. We have pointed out, 
that as a possible side effect of this procedure symmetries 
of the cluster might be destroyed, leading to certain un- 
physical features. 

We intend to use the results of this work as a limit- 
ing case to test various approximations for the complete 
CKLM. Furthermore, the cluster result itself can serve as 
a starting point for an approximation of the lattice. We 
are convinced that techniques similar to what has been 
proposed here can be used to treat slightly more com- 
plex systems, like a cluster with a higher value for the 
localized spin S or another number of lattice sites. 
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APPENDIX A: THE HUBBARD CLUSTER 

An analytical solution of a two-site Hubbard cluster has already been given by several authorsEJBI. We present 
here a direct solution, that is compatible to the language of equations of motion used in the rest of the paper. 

We study the Hamiltonian 



H 



E f( c| 



-2cr L lcr 



2 
k=l 



l k-c 



(Al) 



Using abbreviation (||) we need to consider a set of Green's functions which is almost identical to the list given 
in section y. However, in contrast to the CKLM cluster we have 5=1 what decreases the number of functions 
drastically. We are left with the following set of equations of motion, written again in a matrix representation: 



E 



( Ga\ 

Gb 
G c 

Gd 

Gf 
G g 
G h 
Gj 
Gk 
G L 

\GmJ 



(ef\U 



U eT eT -eT| 

ef -ef ef 
ef -ef ef 
-ef ef ef U 



U 



eT + U 



u 



-u 



U eT eT -eT\ 
ef -ef ef 
ef -ef U ef 
-ef ef ef 



\ 



U 



eT+U J 



( Ga\ 




(k a \ 


Gb 




K b 


G c 




K c 


Gd 




K D 


Gf 




K F 


Gg 




K G 


Gh 




K h 


Gj 




Kj 


Gk 




K K 


G L 




«l 


\GmJ 




\K M ) 



(A2) 
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It is now essential to introduce the combinations 



xf 
4 e) 


= 4 £) - 

= 4 £) - 


f4 £) , 
f4\ 


4 e) 
4 £) 


= 4 £) 
= 4' 


-4 £) 
-4 £) 


4 £) 
4> 


= 4 £) - 
= 4 £) - 


^4 £) < 


4 £) 


= 4 £) 
= 4 £) 


-4 £) 



(A3) 



where "A" stands for "G" or "A". As one can see the matrix representation of the original equations of motion 
contains blocks of the size 4x4. After the proposed combination the remaining maximum block size is 2 x 2, which 
can easily be inverted: 




E 

7=±1 



1 
1 



u 4T 



E-% -jb 



(A4) 



The form of the energy constant b = y (^J + AT 2 demonstrates, that quadratic system sizes are inevitable. We get 
the following matrix for the equations of motion: 



E 



(G A \ 

Gb 

G c 
Gt 
G ( i 
Gg 
Gh 

Gi 
G k 
G\ 
\Gm J 



U~eT2eT 
2ef -ef 



u 

_2_ 



U+eT 



el 



U 

u 



eT + U 





u 




u 




2 




2 


u 




U 




2 


U 
2 


2 


U 
2 


U 




U 




2 




2 





U-eT2eT 
2ef -ef 



U + eT 



cT 



\ 



U 
U 



eT + U J 



( Ga\ 

Gb 
G c 
Gi 
Gd 
Gg 
Gh 

G S 
G k 
G\ 
\Gm / 



(K A \ 

K h 
K c 
K f 
K d 
K G 
K h 

«j 
A k 

K x 

\K M J 



(A5) 



What remains is a straightforward solution of (A5). The expression of G\J can immediately be written down. To 



get G, and Gl in a form with linear energy poles, a single partial fraction expansion is necessary. For G^ an d 



G- one has to apply (A4). And with a few more steps one gets 



G 



(0 




A", 



(<0 



A, 



(e) 



A', 



(0 



K, 



(<0 



A, 



(e) 



A, 



(0 



2(E-eT) 2{E-eT) 2(E 

l + 7¥+7i)4 £) + (l+7¥ ~7i)A, 



(0 



7^)4 £) 



4(£ + eT - 



u 

2 



lb) 



2(E + eT 



u 

2 



• 7 6) 



(A6) 



We have defined the sum of correlation functions 



As a last step on the way to the desired single-particle Green's function one has to undo the combination 
considering the sum | (G^ +) +G { A ] y Then we obtain the result 



by 
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8(E + f-^ + b) 
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APPENDIX B: COMBINATIONS OF GREEN'S FUNCTIONS 



First of all the notation for the Green's functions needs to be clarified. Each of them gets to subindex that consists 
of a number and a letter. The latter part has already been introduced above when the list of Green's functions 
with different possibilities for the fermionic operators is given. The former part enumerates the possibilities for 
spin operators S. The numbers 1, . . . , 10 are used for § = 1, S z s: S§, S'", Sr CT , S%S§, SIS'" , S z s Sg a , Sg a S°, S'" S? . 

Additionally, the Green's functions (^Sc\_ a c Sa c sa j)^ x , and ((«Sn aa Cj -c g _ .c a _ .))p with x = 1,...,5 for S — 
Sg , Sg , SjjS% , S Z S§ , S^ a S~ a appear in the necessary equations of motion. Each Green's function taken 8 times 
for the different possibilities for s, a and e gives the number 1040. 

The first step of combinations is determined by the spin operators. We use the rules that have been applied to the 
insulatonj. 
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H {±) 
n P2 
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The ne xt step produces the correct combinations of fermionic operators of different Green's functions. We use the 
rules ( |A3p applied to the Hubbard cluster. 
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